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Abstract 

We study the scalar condensate and the topological susceptibility for a continuous range of 
quark masses in the Schwinger model with Nf = 0, 1, 2 dynamical flavors, using both the 
overlap and the staggered discretization. At finite lattice spacing the differences between 
the two formulations become rather dramatic near the chiral limit, but they get severely 
reduced, at the coupling considered, after a few smearing steps. 
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1 Introduction 



Recently, the staggered action when coupled to "HYPed" backgrounds [T| has attracted renewed 
^ ■ interest. This is mainly due to the cheapness of this formulation which bears the promise that 
realistic unquenched simulations can be performed with currently available resources j2] . More- 
over, as a remnant of the full SU(Nf = d)A group, the continuous "two-hop global" symmetry 

g. X (x)^exp(i0 A (-l)2>)x(x) , ^)^x(^)exp(-i^(-l) En ") , (1) 

with x = a(rii, . . . , n^) protects the fermion mass against additive renormalization. On the other 
hand, taking the square and quartic root of the determinant (to obtain Nf = 2 + 1 dynamical 
flavors) might spoil the locality of the action [Hj, and the often believed insensitivity to topology 
could undermine attempts to push towards the chiral limit. 

The other extreme in terms of computational cost is represented by the closely related 
domain-wall jl] and overlap [3] fermions. For these actions, the Ginsparg- Wilson relation jB] 
rN ; 

Di 5 + l5 D = 0, 75 = 75(1 D) (2) 

P 

holds in the massless limit (with p for the moment an arbitrary papameter that will be specified 
in ©)> implying invariance under the full chiral symmetry group [7] 

Sip = i 5 ip, Sip = V>75 (3) 

at finite lattice spacing, which again excludes additive mass renormalization and prevents op- 
erators in different chiral multiplets from mixing. 

Given this situation, we decided to investigate the difference between staggered and overlap 
fermions at finite lattice spacing in a simple theory where the concept of chiral symmetry proves 
relevant and some interesting quantities are known analytically. These criteria are fulfilled by 
the generalized Schwinger model (QED in 2D with Nf massive degenerate fermions), a super- 
renormalizable theory jSj where the scale is set through the dimensionful fundamental coupling 
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We produced 10,000 independent gauge configurations on a NxN = 20 2 lattice with the 
standard (compact) Wilson gauge action at (5 = 4, giving a plaquette of 0.86279(10). For 
each configuration we determine the complete eigenvalue spectrum of the massless overlap and 
staggered operators. This allows us to compute - for any given mass - the condensate and the 
fermion determinant, which we use to reweight our observables to Nf = 1,2 [9 . 

We define the massive overlap operator as 

DZ = (1- -tt)D ov + m , aD™ = 1 + 7s sign(a 75J D W p ) (5) 
2p 

with D_ the Wilson operator at negative mass —p/a. We further set p = l, which we checked, 
following Ref. [TJIj, is an almost optimal choice with respect to locality at our coupling. 

Previous work on the Schwinger model using a direct approach for the computation of the 
scalar condensate is collected in [TTl[T2l[TH] for staggered, domain- wall/overlap or both actions, 
respectively. We checked that we reproduce the staggered condensate from [THj and the overlap 
condensate of FHLW/GHR in \T2\. 



2 Scalar condensate 

In the chiral limit, the Nf — 1 scalar condensate is given by [S] 

Xsca(ra = 0) e 7 



g 2tt 3 / 2 



0.1599... [N f = l]. (6) 



For Nf = 2 a. non-zero value would signal spontaneous symmetry breaking and therefore violate 
the Mermin-Wagner-Hohenberg-Coleman theorem . The prediction how explicit symmetry 
breaking modifies this zero is 



= 0.388... (-) 1/3 [N f = 2). (7) 
In the staggered formulation we follow jll] and implement the (1-flavor) condensate through 

Xsca 1 1 



9 " 2N*g 



(XX) , (8) 



where the purpose of the factor 1/2 is to compensate the two-fold degeneracy of the staggered 
formulation in 2D. Denoting the eigenvalues of the massless staggered Dirac operator by A (they 
show up in complex conjugate pairs with zero real part), the reweighted condensate is 



Xsca 1 1 (detTO^E^) 
g 2I?g (det(^)^/ 2 ) 



det(££) = n(A+m) . (9) 



with L = Na and the sum and product running over the entire spectrum. In [TI] one finds also 
the associated free field limit 

o N/2 , 

X E1 = J L _my 1 

g N 2 g ^ (am) 2 + sm(2ni/N) 2 + sm(2nj/N) 2 K ' 

which is, of course, independent of Nf. 

For overlap fermions, the scalar condensate is unambiguously defined as [TJ] 

^ = -^0*' ^ = 1± ?V en) 
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Denoting the eigenvalues of the massless overlap Dirac operator by A, and remembering that 
we work with p= 1, the reweighted condensate is 



— = Ug- (det(D%m ' det ^) = n((l-^-)A+m) (12) 

where the sum runs over the full spectrum. These eigenvalues occur either in complex conjugate 
pairs or as isolated chiral (doubler) modes at aA = (2). Finally, one can rewrite (fT2*J) as 

Y l (detO^)^ E' tM A / a \-i 



g " (det( J D-)7V /) > V 2 

where aA is purely imaginary and the sum excludes the doubler modes at aA = 2. 

Due to the (remnant) chiral symmetry of (staggered) overlap fermions, no subtraction of the 
condensate is required. Furthermore, the Schwinger model is super-renormalizable, i.e. there is 
no renormalization needed at infinite cut-off. In fact, due to the dimensionful coupling (J3J), all 
lattice renormalization factors have the form Z = 1 + 0(a 2 g 2 ) and therefore are 1, up to 0(a 2 ) 
corrections. Hence, both the staggered and the overlap discretization yield (for any Nf and 
m>0) a finite condensate which is subject to 0(a 2 ) cut-off effects, without a need for additive 
or multiplicative renormalization. 

A modification which is motivated by what is done in full QCD simulations, is to consider 
the Dirac operator on a "copy" of an element of the Markov chain on which one applies one or 
more APE/HYP smearing steps (which in 2D is the same). One then thinks of the operator 
as one in the original links, which is less local. For the staggered action, such a modification 
preserves the universality class while considerably reducing the "taste" violation jTHj. Being 
unaware of any detailed understanding of what the optimum smearing parameter is, we decided 
to combine the staple and the original link with equal weight, which means in a £7(1) theory 
that one takes the arithmetic mean of the phases. 

Fig. [I] displays our results for the scalar condensate. As expected, the overlap condensate 
exhibits the qualitatively correct behavior in the chiral limit. For Nf — it shows the quenched 
divergence, while for Nf — 2 it tends to zero indicating the absence of spontaneous chiral 
symmetry breaking. In the Nf = l case it tends to a constant which seems compatible with the 
analytic result (JOJ) in infinite volume, marked by an asterisk. Here, we dare to compare to the 
infinite volume value, since the results in [TJj suggest that for Nf = 1 in a box with L > 1/ (y/ng) 
finite volume corrections are exponentially small. On the other hand, the staggered condensate 
exhibits a qualitatively wrong behavior in the chiral limit, vanishing for all Nf. At large masses 
it differs substantially from the overlap condensate, presumably due to cutoff effects. In the 
intermediate region, there is a smooth turnover with no sign of a quenched divergence. 

The dashed curve is Smilga's infinite volume result (jZJ), while the dash-dotted graph rep- 
resents the free- field expression (jl()[). Finally, since changing m at fixed box- volume V means 
that one moves from the large to the small Leutwyler-Smilga regime, the point where the 
Leutwyler-Smilga parameter [TS] 

x = VT,m (14) 

equals 1 is indicated with a vertical dotted line (S denotes the analytical result (JBJ))- 

The real surprise comes when the operators are evaluated after one step of APE/HYP 
smearing. The overlap condensate stays (for any Nf) virtually unchanged, except that the size 
of the 0(a 2 ) artefacts is modified. Formally, the same statement holds true with the staggered 
action, but from a quantitative point of view the change is rather dramatic. In the Nf = 1 
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Figure 1: Scalar condensate after (top) and 1 (bottom) steps of APE/HYP smearing. 
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Figure 2: Hybrid condensate after (top) and 1 (bottom) steps of APE/HYP smearing. 
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Figure 3: Scalar condensate after (top) and 3 (bottom) steps of APE/HYP smearing 
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case one can "trust" it down to much smaller quark masses, until it finally collapses. And for 
iV} = it shows a nice blow-up at moderately small quark masses. If one really considers the 
limit m — > it still tends to zero, but the point where this happens is deep within the e regime 
(x<Cl) [18J. The effect of additional smearing steps will be discussed below. 

Fig.|21 shows our results for the "hybrid" condensate, labeled after the formulation used 
for the valence quark. The "staggered" condensate (with Nf = 1,2 sea quarks built from the 
overlap determinant) looks virtually unchanged w.r.t. Fig. ^ (the curve with Nf = is identical). 
On the other hand, the "overlap" condensate (with Nf = 1,2 sea quarks constructed from the 
staggered determinant) looks much worse than the original (true) overlap version; there is a 
divergence near m = for any Nf. This shows that the failure of the unsmeared ("naive/thin- 
link") staggered formulation cannot be attributed to either the determinant or the valence 
prescription alone - it's the cancellation of the zero in the determinant and the "one over zero" 
in the propagator which is needed to get the finite value (jBJ) in the 1-flavor case. 

Upon applying just one step of smearing the situation improves dramatically; both "hybrids" 
look qualitatively right down to much smaller quark masses - but still, eventually the non- 
faithful representation of the zero-mode(s) in the staggered part reflects itself in a fake blow-up 
or drop-down close to the chiral limit. 

Fig-IHl contains our results for the condensate in a log-log representation. The overlap data 
show that the analytic structure in a finite volume is indeed 

f 1/m (N f = 0) 
Xscai oc < const (N f = l) (15) 
[ m (N f = 2) 

near the chiral limit and consistent with 

Xscai oc m 1/3 (16) 

at large masses, albeit the coefficient in ((7|) is (presumably due to cut-off effects) not reproduced. 

The bottom part demonstrates that 3 APE/HYP-smearing steps manage to completely 
change the overall picture of the staggered condensate while the overlap condensate gets rescaled 
by a factor which differs from 1 only marginally. Now, the two formulations are in perfect 
agreement at (moderately) small quark masses, showing only a mild discrepancy at larger m/g: 
the staggered condensate seems to tend towards the free-field limit (|10p. while the overlap 
version moves closer to Smilga's prediction (J7J). Eventually, the staggered condensate vanishes 
in the chiral limit - even for Nf = 0. The turnover point, however, is off the scale of the plot. 



3 Selection theorem 

A remarkable feature of the Schwinger model with Nf = 1 is the "selection theorem". This 
theorem states that in the chiral limit the non-zero value JBJ) is formed exclusively from the 
zero-modes of chirality ±1, which live on backgrounds where the topological charge 

q(A) = |- J e^F,!,, dx e Z (17) 

takes the value =pl. This statement alludes to a combination of the Atiyah-Singer index theorem 
and the vanishing theorem. The former relates ()17|) to the index of the Dirac operator, defined 
as the difference of the number of negative and positive chirality zero-modes 

ind(A) = n_ - n + , (18) 
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in a simple manner 

q(A) = md(A) , (19) 

while the latter says that in 2D there is no configuration which supports both positive and 
negative chirality zero-modes [THj 

n + = , . 

n+ ^ n_ = ' { > 

To prove the "selection theorem" one starts from the partition function of the massless theory 
in the finite volume with periodic boundary conditions (Nf>l) J7| 

/ |,| _ \ Nf 

Z[fj,r]) = NY; IDAtoe-M** m^Xfe) det ' (D) N f e + 1 ^ (21) 

gGZ J \k=l J 

from which the condensate is computed by taking one derivative w.r.t. r\ and r\ and then setting 
the external sources to zero. The outer sum is over all topological sectors and the product is 
over the \q\ positive or negative chirality zero-modes. The primes indicate that the determinant 
and the fermion Green's function are computed on the subspace orthogonal to the zero-modes. 
A non-zero value is obtained only if both derivatives hit the prefactor, leaving nothing but the 
explicit zero-mode and its conjugate behind. This can happen only for Nf = l and results in a 
condensate which is generated exclusively by the zero-modes with chirality ±1. Note that the 
argument will hold true in QCD, too, once the current body of numerical evidence in favor of 
a vanishing theorem in 4D (see Fig. 1 in [20]) has been replaced by a mathematical proof. 

As a technical point we mention that on the lattice we define the topological charge of a 
background U as the index of the massless overlap operator [2U Ej 

q(U) = ind(f/) = ^tr(a 7 5D ov ) (22) 
and use it for both the staggered and the overlap evaluation of the sectoral condensate. 

Fig. |3] displays our results for the condensate if we truncate the partition function at a given 
1 9 i max- The overlap construct faithfully reproduces the "selection theorem", which means that 
the Nf — 1 condensate tends to zero if one restricts it to the topologically trivial sector, while 
it takes the Schwinger value (J0J) in the chiral limit for any other |g| max . 

Analogously, the staggered condensate exhibits (for both Nf shown) a many-sigma difference 
between |g| max = 0, 1,2. This disproves the widely believed fiction that staggered fermions are 
"insensitive to topology" - but they are not sensitive the right way; they don't seem to know 
about the index and the vanishing theorems which are at the root of the "selection theorem" . 

Fig.EJshows that the real surprise comes again after one or a few APE/HYP-smearing steps. 
Already one step lets the staggered Nf = l condensate (at intermediate mass) develop a marked 
sensitivity on the topological charge of the background. 

After two more steps the qualitative picture is just like in the (unsmeared) overlap case (c.f. 
Fig.HJ), i.e. smeared staggered fermions do know about the relationship between topology and the 
chirality of zero-modes - down to rather small (but non-zero) quark masses. 

Of course, if one really performs the chiral limit, the staggered condensate still tends to zero 
- for any |g| max and even in the quenched case (Nf = 0). This is visible in the top of Fig.0 
while in the bottom part it is off the scale. 

Since in the Schwinger model various observables have been seen to depend on the topolo- 
gical charge j22], it is surprising that the "selection theorem" has not been checked before. 
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Figure 4: Sectoral condensate for overlap (top) and staggered quarks without smearing. 
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Figure 5: Sectoral condensate for staggered quarks with 1 and 3 steps of APE/HYP smearing. 
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4 Topological susceptibility 



Another interesting observable to study the effects of dynamical fermions is the topological 
susceptibility which, in the context of this note, shall be defined (in the continuum) through 

(det(P+m)"/ q 2 ) 
X ^ = ^ V (det(D +m rf) ■ (23) 

The main difference to the scalar condensate is that the topological susceptibility depends only 
on the sea-quarks, thus offering a potentially cleaner view at the effects of square-rooting the 
staggered determinant to get Nf = l. 

For staggered quarks, the definition (j23J, taken in fixed volume, reduces to 

Xtop l_ (det(D%) N f/ 2 q 2 ) 

g 2 N 2 (det(D%) N f/ 2 ) ' 1 ' 

while for overlap quarks, the implementation reads 

Xtop _ P (det(DST)^ q 2 ) 



g 2 N 2 (det{D%) N f) 



(25) 



where det(-D^) and det(D™) are defined in H2jl. In either case the sum is over the spectrum 
of the massless operator, and we apply the overlap definition of the topological charge 
both in (j21|) and in 

Finally we like to mention the continuum prediction how the topological susceptibility in 
QCD (!) tends to zero, if the quark mass does [IB] (the LS parameter x was defined in (JTlj) ). 

Xtop = (iV / = lVx»l) (26) 

Xtop oc m N f (x<l). (27) 



Fig.El contains our results for the mass dependence of the topological susceptibility, defined 
via (J23| l2*K|) . The full and the dashed lines represent the prediction of 1)26)1 for the cases 
Nf = 1,2, respectively One can see that for both discretizations Xtop gets suppressed by 
dynamical fermion effects, but close to the chiral limit only the overlap determinant leads to 
results which are compatible with the QCD prediction (J26)) for Nf — 1 and ()27)) for Nf = 2. 

Again, just one smearing step proves sufficient to almost eliminate the lattice artefacts the 
"unsmeared/naive" staggered determinant was plagued with, while the topological susceptibil- 
ity for overlap fermions stays basically invariant under such a modification. 



5 Spectral hint 

As a hint of what is the likely reason behind the remarkable success of one or a few APE/HYP 
smearing steps at the /?-value considered, we like to present the effect of such a modification on 
the spectrum of some individual configurations. In Fig.Qwe plot the physically relevant part 
of the spectrum of the staggered (the a\ from 0) and overlap (the aX from (JEJ)) operator on 
four selected configurations before and after smearing. 

The two configurations on the left are typical examples for topological charge q = and 
\q\ = 1, respectively. One sees that on the unsmeared configurations the staggered spectrum 
does not resemble that of the overlap operator. Furthermore, it is hard to see a qualitative 
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Figure 6: Topological susceptibility with (top) and 1 steps of APE/HYP smearing. 
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difference between q = and |?| = 1 in the spectrum of the original staggered operator. After 
a few smearing steps this picture changes. The eigenvalues of the staggered operator seem to 
form near-degenerate pairs which sit close to a single overlap aX on the same configuration. In 
particular, in the \q\ — 1 case a pair of eigenmodes moves very close to the real axis. Clearly, such 
a shift mimics the effect of the true zero-mode in the overlap counterpart down to rather small 
quark masses. It is only when the mass becomes comparable to the smallest (smeared staggered) 
eigenvalue that the absence of an exact zero-mode matters and a qualitative difference between 
the staggered and overlap Dirac operators (on that configuration) shows up. 

In the third column of Fig.Q a typical configuration with higher topological charge (\q\ =4) 
is displayed. Here, the picture is not so nice any more, since after 3 smearing steps only 3 pairs 
of eigenmodes have come close to the real axis and the fourth one is still somewhat further out. 

To show that these findings are not generic, the last column of Fig.[7| presents a selected 
configuration on which the topological charge varies repeatedly under subsequent smearing 
steps. Qualitative resemblance between the spectra of the two operators is vague at best and 
there is no evidence of a pair of staggered eigenvalues moving close to the real axis. Obviously, 
such "sick" configurations will occur more frequently at larger coupling. 



6 Summary 

Our findings may be summarized by the following statements: 

1. At finite spacing and in a finite box, the "naive" staggered action leads to a scalar con- 
densate which is qualitatively wrong: the staggered results vanish in the chiral limit for 
any Nf, while the overlap successfully reproduces the quenched ~ 1/m divergence, and 
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the analytically known Schwinger value in the chiral limit for Nf = l. 

2. Considering both types of "hybrid" formulations (staggered valence quarks with overlap 
sea quarks and vice versa) we find that the failure of the naive staggered formulation 
cannot be attributed to either the determinant or the propagator alone. 

3. The "selection theorem" is reproduced, in an impressive manner, with overlap fermions, 
while naive staggered fermions fail completely. 

4. The topological susceptibility shows lattice artefacts which are large for unimproved stag- 
gered fermions, while the overlap results seem consistent with the known chiral behavior. 

5. Taking the square root of the staggered determinant to have Nf — 1 seems to be no more 
or less harmful than dividing by 2, in the valence sector, to get the 1-flavor condensate. 
All naive staggered results seem to only gradually vary for Nf = 0, 1, 2 - even near the 
chiral limit, where they shouldn't. 

6. At the f3 value considered, already one smearing step brings a remarkable improvement: 
although formally the chiral limit of staggered fermions is still wrong, the mass at which 
one begins to see this is dramatically lowered. At moderate quark masses, one observes 
a clear blow-up of the APE/HYP smeared staggered condensate for Nf = 0, while it stays 
close to the Schwinger value © for Nf = l and keeps the qualitatively correct behavior at 
Nf = 2. Moreover, in spite of naive staggered fermions being misguided in the way they 
see topology, the APE/HYP variety condensate shows a remarkable sensitivity on |g| max 
for Nf = and Nf — 1, reproducing the "selection theorem" in the latter case. Finally, the 
staggered artefacts in the topological susceptibility get drastically reduced, resulting in 
good agreement with the overlap curve down to very small quark masses. 

7. The main difference after a few smearing steps is the size of lattice artefacts in the 
condensate at large quark masses which make the smeared overlap condensate (for any Nf) 
lie closer to the analytical prediction (|7J) by Smilga, while the smeared staggered values 
move more towards the free field prediction (fTU|) . The qualitative failure of staggered 
fermions in the chiral limit is only visible at very small quark masses. 

8. Under moderate smearing the modes of the massless staggered Dirac operator form near- 
degenerate pairs which move close to the corresponding overlap (single) eigenvalues (A), 
supporting the square root and factor 1/2 procedure for staggered spectroscopy. On non- 
trivial backgrounds the mismatch between the staggered fake zero-mode and the true 
overlap zero gives an estimate of the quark mass at which the deficiency of the staggered 
formulation gets visible, though there are configurations which break this analogy. 

As a warning against an overly optimistic interpretation of 6/8, we feel obliged to recall that 
we were working with fairly smooth gauge fields. It is not clear whether such a nice pattern 
holds true also in QCD at accessible couplings. 
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